TWO CORES OF A NONNEGATIVE MATRIX 
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Abstract. We prove that the sequence of eigencones (i.e., cones of nonnegative eigen- 
vectors) of matrix powers is periodic both in max algebra and in nonnegative linear 
algebra. Using an argument of Pullman, we also show that the Minkowski sum of the 
eigencones of A* is equal to the core of A defined as the intersection of nonnegative 
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column spans of matrix powers, also in max algebra. Based on this, we describe the set 
of extremal rays of the core. 

The spectral theory of matrix powers and the theory of matrix core is developed in 
max algebra and in nonnegative linear algebra simultaneously wherever possible, in order 
to unify and compare both versions of the same theory. 
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1. Introduction 



00 " The nonnegative reals M + under the usual multiplication give rise to two semirings with 

O 

addition defined in two ways: first with the usual addition, and second where the role of 
addition is played by maximum. We refer to the first semiring as nonnegative algebra, 
^ ■ and to the second one as max algebra. Thus we consider the properties of nonnegative 
matrices with entries in two semirings, the semiring of nonnegative numbers with usual 
addition and multiplication, and the semiring called "max algebra" . 

Our chief object of study is the core of a nonnegative matrix A. This concept was 
introduced by Pullman in [33j, and is defined as the intersection of the cones spanned 
(generated) by the columns of matrix powers A 1 . Pullman provided a geometric approach 
to the Perron-Frobenius theory of a nonnegative matrix A based on the properties of the 
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core of A. He investigated the action of a matrix on its core showing that it is bijective 
and that the extremal rays of the core can be partitioned into periodic orbits. In other 
words, extremal rays of the core of A are nonnegative eigenvectors of the powers of A 
(associated with nonnegative eigenvalues). 

One of the main purposes of the present paper is to extend Pullman's core to max 
algebra, thereby investigating the periodic sequence of eigencones of max-algebraic matrix 
powers. However, following the line of [7J (TTJ [23], we develop the theory both in max 
algebra and in nonnegative algebra, in order to emphasize common features as well as 
differences, to provide general (simultaneous) proofs where this is possible, and even to add 
a bit of further information in the case of nonnegative algebra. Our motivation is closely 
related to the Litvinov-Maslov correspondence principle [25], viewing the idempotent 
mathematics (in particular, max algebra) as a "classical shadow" of the "traditional" 
mathematics over real and complex fields. 

Pullman's core can be also seen as closely related to the limits of powers of nonnegative 
matrices. However it is a different concept. Consider the simple example 



Then, for any nonnegative x, A k x will tend to a multiple of (1, 0) T while the core of A is 
the entire nonnegative orthant R+. 

To the authors' knowledge, the core of a nonnegative matrix has not received much 
attention. Tam and Schneider [13] extended the concept of core to cone-preserving maps. 
On the other hand, there is much more literature on the related but distinct question 
of the limiting sets of homogeneous and non-homogeneous Markov chains in nonnegative 
algebra; see the books by Hartfiel [21] and Seneta [37] and, e.g., the works of Chi p2] 
and Sierksma [12] . In max algebra, see the results on the ultimate column span of matrix 
powers for irreducible matrices [6] Theorem 8.3.11 or [38], and by Merlet [27] on the 
invariant max cone of non-homogeneous matrix products. 

The theory of the core relies on the behaviour of matrix powers. In the nonnegative 
algebra, recall the works of Friedland-Schneider [TH] and Rothblum- Whittle [34J (the role 
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of distinguished classes which we call "spectral classes" , algebraic and geometric growth 
rates, applications). The theory of max- algebraic matrix powers is similar. However, 
the max-algebraic powers have a well-defined periodic ultimate behaviour starting after 
sufficiently large time. This ultimate behaviour has been known since the work of Cohen et 
al. [H] (irreducible case), and is described in greater generality and detail, e.g., by Akian, 
Gaubert and Walsh [T], Gavalec [2D], De Schutter [SB], and the authors P ESI BE] of the 
present paper. In particular, the Cyclicity Theorem of Cohen et al. (see [2j El [13], 122] ) 
implies that extremals of the core split into periodic orbits for any irreducible matrix (see 
Subsection 4.2). 

Let us remark that some results on the eigenvectors of max-algebraic matrix powers 
have been obtained by Butkovic and Cuninghame-Green [BJ IE]- The present paper also 
aims to extend and complete the research initiated in that work. 

This paper is organized as follows. In Section |2] we introduce the basics of irreducible 
and reducible Perron-Frobenius theory in max algebra and in nonnegative linear algebra. 
In particular see the notations table in Subsection 12.61 Necessary elements of the graph 
exponents (Boolean cyclicity) theory and the max-algebraic Cyclicity Theorem are also 
given here in Section [2] 

In Subsection 13 . 1 1 we explain the relation between spectral classes of matrix powers, and 
how the eigencones associated with general eigenvalues can be reduced to the case of the 
greatest eigenvalue, see in particular Theorems 13.51 and 13.61 This subsection is very close 
to what is explained, e.g., by Frobenius [17] and, more recently, by Tarn-Schneider [33] 
Section 4 (in the case of nonnegative algebra). In Subsection 13. 2l we show that the sequence 
of eigencones of matrix powers in max algebra is periodic, and compare this result with 
the case of nonnegative matrix algebra, see Theorem 13.91 for the main formulation. Note 
that the nature of this periodicity appears to be different, with no simultaneous proof. In 
Subsection 13.31 the inclusion relation on eigencones is studied further. Namely, we verify 
that these relations precisely follow divisibility, see Theorem 13.121 

In Section H] we show that the core of a matrix A in both algebras equals the sum of the 
eigencones of its powers A k , which equals the (Minkowski) sum of eigencones of a certain 
matrix power (the exponent of this power is different in max algebra and in nonnegative 
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linear algebra). Note that the main argument of Subsection 14.31 only slightly generalizes 
that of Pullman [33], Section 2, to work in both algebras simultaneously. However, the 
main Theorems 14.61 and 14.71 of Section H] describing the matrix core and its extremal 
rays in both algebras (with nonnegative part referring to jl3], Section 4) are based also 
on the periodicity of eigencone sequences shown in Section [31 with essential difference 
between the two algebras. To this end, in the case of max algebra there are wide classes 
of matrices, where the above mentioned result follows from the ultimate periodicity, see 
Subsection 14.21 

Theorem 14.61 is illustrated by a pair of examples in Subsection 14.41 In the first exam- 
ple, the sequence of max-algebraic column spans converges to the (Minkowski) sum of 
eigencones of a certain power only in the limit (just as in the usual nonnegative algebra). 
However the second example shows that in some special cases, finite stabilization of the 
core may take place almost immediately. 

2. Spectral theory: max-algebraic and nonnegative 

2.1. Nonnegative algebra. In this paper we are concerned only with nonnegative eigen- 
values and nonnegative eigenvectors of a nonnegative matrix. In order to bring our termi- 
nology into line with the corresponding theory for max algebra we use the terms eigenvalue 
and eigenvector in a restrictive fashion appropriate to our semiring point of view. Thus 
we shall call p an eigenvalue of a nonnegative matrix A (only) if there is a nonnegative 
eigenvector x of A for p. Further x will be called an eigenvector (only) if it is nonnegative. 

In this paper we are led to state the familiar Perron-Frobenius theorem in slightly 
unusual terms: An irreducible nonnegative matrix A has a unique eigenvalue denoted 
by p + (A), which is positive (unless A is the lxl matrix 0). Further, the eigenvector 
x associated with p + (A) is essentially unique, that is all eigenvectors are multiples of x. 
The nonnegative multiples of x constitute the cone of eigenvectors (in the above sense) 
V + (A,p + (A)) associated with p + (A). 

A general (reducible) matrix A e M+ Xra may have several positive eigenvalues with asso- 
ciated cones of nonnegative eigenvectors (eigencones) , and p + {A) will denote the biggest 
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such eigenvalue, in general. Eigenvalue p + (A) is also called the principal eigenvalue, and 
V + (A, p + (A)) is called the principal eigencone. 

Recall that a subset V C R™ is called a (convex) cone if 1) av G V for all v G V and 
a G R+, 2) w + t> G V for u,v £ V . Note that cones in the nonnegative orthant can be 
considered as "subspaces" , with respect to the semiring of nonnegative numbers. In this 
vein, a cone V is said to be generated by S C R™ if each i> G V can be represented as a 
nonnegative combination v = XLeS a ^ x where only finitely many a x G R+ are different 
from zero. A vector z in a cone V is called an extremal, ii z — u + v and ti,«6 V imply 
z = = a v v for some scalars a u and cv Any closed cone in R™ is generated by its 
extremals; in particular, this holds for any finitely generated cone. 

Let us recall some basic notions related to (ir)reducibility, which we use also in max 
algebra. With a matrix A = (a^) G R™ xra we associate a weighted (di)graph Q(A) with 
the set of nodes N = {1, . . . ,n} and set of edges E <0 N x N containing a pair (i,j) if 
and only if ^ 0; the weight of an edge G E is defined to be w(i,j) := a^. A 
graph with just one node and no edges will be called trivial. 

A path P in Q(A) consisting of the edges (i , i±), (ii, i 2 ), ■ ■ ■ , (h~i,h) has length I (P) := t 
and weight w(P) := w(i , ii) ■ w(ii, i 2 ) • • -w{i t -\, it), and is called an i — j path if i — i 
and i t = j. P is called a cycle if io — it- 

Recall that A = (a^) G IR™ xn is irreducible if G(A) is trivial or for any i, j G {1, . . . , n} 
there is an i — j path. Otherwise A is reducible. 

Notation A xfc will stand for the usual fcth power of a nonnegative matrix. 

2.2. Max algebra. By max algebra we understand the set of nonnegative numbers M + 
where the role of addition is played by taking maximum of two numbers: a@b := max(a, b), 
and the multiplication is as in the usual arithmetics. This is carried over to matrices and 
vectors like in the usual linear algebra so that for two matrices A = (ay) and B = (by) 
of appropriate sizes, (A © B) v - = © b^ and (A © B) ik = @ k a ik b k j. Notation A® k will 
stand for the /cth max-algebraic power. 

Max-linear spaces are analogues of convex cones. A set V C will be called a max 
cone if 1) av G V for all v G V and a G R+, 2) u@ v G V for u, v G V. Max cones are also 
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known as idempotent semimodules, see (251 [22]- A max cone V is said to be generated 
by S C if each v E V can be represented as a max combination t> = ® xeS ct x x where 
only finitely many (nonnegative) a x are different from zero. When V is generated (we 
also say "spanned") by S, this is denoted V = span e (S'). When V is generated by the 
columns of a matrix A, this is denoted V = span e (A). 

A vector z in a max cone V C IR™ is called an extremal ii z — u @ v and u,v E V 
imply 2; = ti or 2 = w. Any finitely generated max cone is generated by its extremals, see 
Wagneur [H] and [H [19] for recent extensions. 

The maximum cycle geometric mean of A is defined by 

(1) X(A) = max{w(C) W) ; C is a cycle in 0(A)} 

(see Sect. 12.11 for definitions of Q(A) and cycles). The critical graph of A, denoted by 
C(A), consists of all nodes and edges belonging to the cycles which attain the maximum 
in ([!]). The set of such nodes will be called critical and denoted N c ; the set of such edges 
will be called critical and denoted E c . Observe that the critical graph, defined as above, 
consists of several strongly connected subgraphs of G(A). Maximal such subgraphs are 
the strongly connected components of C(A). 

If for A G M" xn we have A (g) x = px with p 6 M+ and a nonzero x G , then p is a 
max(- algebraic) eigenvalue and a; is a max (-algebraic) eigenvector associated with p. The 
set of max eigenvectors x associated with p, with the zero vector adjoined to it, is a max 
cone further denoted by Vq(A,p). 

An irreducible A G M" xn has a unique max-algebraic eigenvalue equal to the m.c.g.m. 
A (A) [21 El HH [22]- In general A may have several max eigenvalues, which we will describe 
later. The greatest such eigenvalue will be denoted by p®(A), and called the principal 
max eigenvalue of A. In the irreducible case, the unique max eigenvalue p®(A) = X(A) 
will also be called the max(- algebraic) Perron root. When max algebra and nonnegative 
algebra are considered simultaneously (e.g., Sect. 13.11 or fl~3l) . the principal eigenvalue is 
denoted by p(A). 
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Unlike in nonnegative algebra, there is an explicit description of V©(A, p®(A)), see 
Theorem 12 . 1 1 1 b elow . This description uses the Kleene star 

(2) A* = 1 © A © A m © A m © ... , 

where I denotes the identity matrix. Series fl2]) converges if and only if p®(A) < 1, in 
which case A* = I © A © ... © A^ n ~^ . Note that if p®(A) ^ 0, then p®(A/p®(A)) = 1, 
and (A/ p®(A))* always converges. 

The path interpretation of max-algebraic matrix powers A® 1 is that each entry afj is 
equal to the greatest weight of i — j paths with length I. Consequently, for i ^ j, the 
entry a*- of A* is equal to the greatest weight ofi — j paths (with no length restrictions) . 

2.3. Cyclicity and visualization in max algebra. For a strongly connected graph Q 
define its cyclicity a as the gcd of the lengths of all elementary cycles and the cyclicity 
of a trivial graph to be 1. It is known that if there is an % — j path, then there exists a 
number I such that l(P) = /(mod a) for all i — j paths P. 

When the length of an % — j path is a multiple of a (and hence we have the same for 
all j — i paths), i and j are said to belong to the same cyclic class. When the length of 
this path is 1 modulo a (in other words, if l(P) — 1 is a multiple of a), the cyclic class of i 
(resp., of j) is previous (resp., next) w.r.t. the class of j (resp., of i). The cyclic classes of 
Q form a simple cycle of length a, in terms of these adjacency relations. See [6], Chapter 
8 and jH [381 EH] f° r more information. Cyclic classes are also known as components of 
imprimitivity jl]. 

We define the cyclicity of a (general) graph containing several maximal strongly con- 
nected components (such as the critical graph C(A)) to be the 1cm of the cyclicities of the 
components. For a graph Q = (N, E) with N — {1, . . . , n}, define the associated matrix 
A = (a.y) G {0, l} nxn by = 1 <^ (i, j) G E. This is a matrix over the Boolean semiring 
B := {0, 1}, where addition is the disjunction and multiplication is the conjunction op- 
eration. This semiring is a subsemiring of max algebra, so that it is possible to consider 
the associated matrix as a matrix in max algebra whose entries are either or 1. 

If a Boolean matrix A is associated with a graph Q, then by Q k we denote the graph 
associated with the Boolean power A k . We will need the following observations. Powers 
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of Boolean matrices (over the Boolean semiring) are a topic of independent interest, 
see Brualdi-Ryser [4], Kim [24]. Recall that Q = (N,E) is complete if {i,j) G E for all 
i,j G N. 

Theorem 2.1 (Boolean Cyclicity [24]). Let Q = (JV, E), with the set of nodes N = 
{1, . . . , n}, be strongly connected with cyclicity o. Then 

(i) there exists T(Q) < (n - l) 2 + 1 such that g t+a = Q l for all t > T(G), 

(ii) for t > T(Q) and t a multiple of a, Q l consists of a complete subgraphs not 
accessing each other. 

Inspired by the combinatorial matrix theory [4], the threshold T(Q) can also be called 
the graph exponent. 

Theorem 2.2 (cf. [4J Theorem 3.4.5). Let Q be a strongly connected graph with the set 
of nodes N and cyclicity a. 

(i) Q k consists of gcd (k, a) strongly connected components not accessing each other. 
The union of their sets of nodes is N. 

(ii) Each of the components contains a/(gcd(k,a)) cyclic classes of Q . 

Corollary 2.3. Let Q be a strongly connected graph with cyclicity a, let k, I > 1, then 
gcd(k, a) divides gcd(l, a) if and only if Q k and Q l are such that the node set of every 
component of Q l is contained in the node set of a component of Q k . 

Proof. "If". All components of Q k contain an equal number of cyclic classes, and these 
cyclic classes are distributed among the components of Q . All components of Q l contain 
an equal number of cyclic classes, implying that each component of Q k contains the same 
number of components of Q l . This shows that gcd(/c, a) divides gcd(/,cr). 

"Only if." Observe that the node sets of Q k and Q^ k ^ (or Q l and £s cd C> CT )) are the 
same: since gcd (A;, a) divides k, each component of Qscd(k,a) S p\[^ s j n ^ seV eral components 
of Q k , but the total number of components is the same, hence their node sets are the 
same. The claim follows since each component of £g cd ( fc >°") splits into several components 

of £gcd(i,a) □ 
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There are two max-algebraic extensions of Theorem 12.11 

Theorem 2.4 (Cyclicity Theorem, Cohen et al. [13J). Let A G M™ xn be irreducible and 
let a be the cyclicity ofC(A). Then there exists T(A) such that A® {t+CT *> = {p®(A)) a A® t 
for allt> T(A). 

Theorem 2.5 (Cyclicity of Critical Part, Nachtigall [30]). Let A G M+ Xn and let a be 

the cyclicity ofC(A). Then there exists T C (A) < n 2 such that Af t+a) = (p®(A)) a Af and 
A ®(t+o) = ( p( B( A ^* A ft for all • e N ^ and f > T ^ A y 

Theorem 12.41 is standard [2l[6l[22], and Theorem 12.51 can also be found as [6], Theorem 
8.3.6. Here A^., resp. Aj denote the zth row, resp. the ith column of A. 

Cyclicity in max algebra is related to the Boolean cyclicity by means of the diagonal 
similarity scaling. If x G W] is positive, i.e., it has all components nonzero, then we can 
define X = diag(x) as the diagonal matrix whose diagonal entries are equal to the cor- 
responding entries of x, and whose off-diagonal entries are zero. The diagonal similarity 
scaling A i— > X~ l AX does not change the weights of cycles and eigenvalues (both nonneg- 
ative and max); if z is an eigenvector of X -1 AX then Xz is an eigenvector of A with the 
same eigenvalue. This scaling does not change the critical graph C(A) = (N C (A) , E C (A)) . 
Observe that (X~ x AX) k = X~ 1 A k X (both max and nonnegative) , also showing that 
the cyclicity thresholds T(A) and T C (A) do not change after scaling. The technique of 
nonnegative scaling can be traced back to the works of Fiedler-Ptak |15j . 

Theorem 2.6 (Strict Visualization PI4T]). For each A e M+ Xn with m.c.g.m. \{A) = I, 
there exists x e such that for X := diag(x) ; the matrix B = X~ X AX has b{j < 1, and 
moreover, bij = 1 if and only if G E C (A)(= E C (B)). 

Scaling of Theorem l2.6l is called the strict visualization scaling, and matrix B satisfying 
b^ < 1 and b^ = 1 G E C (B) is called strictly visualized. 

If A = (ay) has all entries < 1, then it makes sense to consider the Boolean matrix 
A^ with entries 
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If A has all entries m, < 1 then = (A^)® k . Note that if A is strictly visualized 

then AM is the associated matrix of C(A). In particular, it follows that any power of a 
strictly visualized matrix is also strictly visualized, and that C{A® k ) = (C(A)) h . After 
this property is observed for strictly visualized matrices, it is extended to the general case 
by means of the strict visualization scaling (Theorem 12 .6 j) . Using C(A® k ) = (C(A)) k , we 
see that T C (A) > T(C(A)), where T{C(A)) is the greatest graph exponent of the strongly 
connected components of C(A), and also that a is the least (i.e., exact) ultimate period 
in Theorems 12.41 and 12.51 

Similarly if a vector x G has Xi < 1, we define x^ having xf = 1 if Xi = 1 
and xf^ = otherwise. Obviously if A and x have all entries not exceeding 1 then 
(Ac)M = A^xM. Naturally supp(xM) is the set of all indices i where x]p = 1. 

Let us conclude this section with a well-known number-theoretic result concerning the 
coin problem of Frobenius, which we see as basic for both Boolean and max-algebraic 
cyclicity. 

Lemma 2.7 (e.g.,jl] Lemma 3.4.2). Let rii, . . . ,n m be integers such that 

gcd(rii, . . . , n m ) = k. Then there exists a number T such that for all integers I with kl > T , 

we have kl = tini + . . . + t m n m for some t 1; . . . , t m > 0. 

2.4. Frobenius normal form. Every matrix A = (aij) G M™ xn can be transformed by 
simultaneous permutations of the rows and columns in almost linear time to a Frobenius 
Normal Form (FNF) P g] 

... \ 
A 22 ... 
... A^ ... 

A r 2 ... A rr J 

where An, A rr are irreducible square submatrices of A, corresponding to the partition 
A^i U . . . U N r — N. Generally, Akl denotes the submatrix of A extracted from rows with 
indices in K C iV and columns with indices in L C N, and A^ is a shortcut for An n v . 
Accordingly, the subvector xjv m of x with indices in iV M will be written as x M . 



(4) 



(A 



ii 

-421 



V A rl 
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Consider the graph associated with A^. It is the same as the graph induced by N^. 
Q{A m ) = (Np, E R (Np x Nfj)). It follows that each of the graphs G(A^) (ji = 1, r) is 
strongly connected and an arc from iV M to N v in Q(A) may exist only if p > v. 

If A is in the Frobenius Normal Form (J4]) then the reduced graph, denoted R(A), is 
the (di)graph whose nodes correspond to N^, for p = l,...,r, and the set of arcs is 
{(//, u);(3keN ll )(3£eN v )a k t>0}). 

In max algebra and in nonnegative algebra, the nodes of R(A) are marked by the 
corresponding eigenvalues (Perron roots), denoted by p® := p®(A MM ) (max algebra), 
p® := p + (A flfl ) (nonnegative algebra), and by p M when both algebras are considered 
simultaneously These nodes will be called classes of A. We naturally attribute to a class 
p also the graph Q(A^) with set of nodes iV M and cyclicity <r M . We will also sometimes 
refer to as a class of A, by abuse of terminology 

We say that a class p is trivial if A w consists of a single diagonal zero entry, i.e., when 
5(j4 w ) is trivial. Class p accesses class v, denoted p — > v, if there exists a p — z/ path 
in i?(v4). A class is called initial, resp. /maZ, if it is not accessed by, resp. if it does not 
access, any other class. Node i accesses class u, denoted by i — > v, if i belongs to a class 
H such that p — > v. 

Simultaneous permutations of the rows and columns of A are equivalent to calculating 
P~ X AP, where P is a permutation matrix. Such transformations do not change the 
eigenvalues, and the eigenvectors before and after such a transformation only differ by 
the order of their components. Hence we will assume without loss of generality that A is 
in FNF dU). 



2.5. Reducible spectral theory. In this section we show two versions of the Perron- 
Frobenius reducible spectral theory. Almost all results are standard: the nonnegative 
part goes back to Frobenius [T7] . Sect. 11, and the max-algebraic counterpart is due to 
Gaubert [H], Ch. IV (also see [6 J for other references). 
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A class v of A is called a spectral class of A associated with eigenvalue p ^ 0, or 
sometimes (A, p)-spectral class for short, if 

p® = p, and p — > v implies p® < p® (max algebra), 

(5) 

p + = p, and p, — > v implies p® < p® (nonnegative algebra). 

In both algebras, note that there may be several spectral classes associated with the 
same eigenvalue. In nonnegative algebra, spectral classes are rather called "distinguished 
classes" [35J, taking into consideration the generalized eigenvectors. We do not consider 
generalized eigenvectors here (and no suitable max-algebraic analogue of this notion is 
known to us). 

Denote by A + (A), resp. A e (A), the set of nonzero eigenvalues of A e M™ xn in non- 
negative linear algebra, resp. in max algebra. This will be denoted by A(A) when both 
algebras are considered simultaneously. The following description is standard. 

Theorem 2.8 Q6J Th. 4.5.4, [35] Th. 3.7). Let A e M+ Xn . Then 

1. A + (A) = {p® 7^ 0; v is spectral} = {p® ^ 0; V/i p, — > v implies p+ < p+} 

2. A e (A) = {p® j£ 0; v is spectral} = {p® ^ 0; Vp p v implies p® < p®} 

In nonnegative algebra, [35], Th. 3.7 describes both spectral classes and eigencones 
associated with general eigenvalues immediately. Proposition 12.91 and Theorem 12.101 can 
also be seen as its corollaries in the case of nonnegative algebra, see below. However, we 
prefer to split the formulation, following the exposition of [6] . An alternative simultaneous 
exposition of spectral theory in both algebras can be found in [23] . 

In both algebras, for each p e A (A) define 

-i f° \ 

A p := p \ , where 

(6) ' \0 A MpMp ) 

M p := {i; i — > u, v is (A, p)-spectral} . 

The next proposition, holding both in max algebra and in nonnegative algebra, allows 
us to reduce the case of arbitrary eigencone to the case of principal eigencone. We will 
assume an appropriate ordering of indices (i.e., nodes of the corresponding graph). Such 
assumptions will sometimes be made also in the sequel, without special mention of them. 
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Proposition 2.9 (|Sl HSl EE])- For A G M+ Xn and each p G A{A), we have V{A,p) = 
V(A p , 1), where 1 zs the principal eigenvalue of A p . 

In max algebra, using Proposition 12. 9[ we define the critical graph associated with 
p G A® (A) as the critical graph of Am p - The strongly connected components of critical 
graphs associated with p, for all p G A e (A), will be called the critical components of A. 
We will see that, in max algebra, the role of spectral classes of A is rather played by 
these critical components, which will be (in analogy with classes of FNF) denoted by p, 
with the node set Np_. In analogy with "//" , we will later use the shortcuts A^ for the 
submatrix A n at., and x^ for the subvector xn~. 

We now describe the principal eigencones in nonnegative linear algebra and then in 
max algebra. By means of Proposition 12.9} this description can be obviously extended to 
the general case. Here, the full-size vectors x^ and x^ are not to be confused with the 
subvectors x M and x^ introduced above. 

Theorem 2.10 (Frobenius- Victory [35] Th. 3.7). Let A G M+ Xn have p + (A) = 1. 

(i) Each spectral class p with p+ = 1 corresponds to an eigenvector x^, whose support 
consists of all indices in all classes that have access to p, and all vectors x of 
V + (A, 1) with suppa; = suppx^ are multiples ofx^K 

(ii) V + (A,1) is generated by x^> of (i), for p ranging over all spectral classes with 

(hi) a;W of (i) are extremals ofV+(A, 1). (Moreover, x^ are linearly independent.) 

Note that the extremality and the usual linear independence of (involving linear 
combinations with possibly negative coefficients) can be deduced from the description 
of supports in part (i), and from the fact that in nonnegative algebra, spectral classes 
associated with the same p do not access each other. This linear independence also means 
that V+(A, 1) is a simplicial cone. See also [35] Th. 4.1. 

Theorem 2.11 (© Th. 4.3.5, [2] Th. 2.8). Let A G M^ xn have p®{A) = 1. 

(i) Each component p ofC(A) corresponds to an eigenvector x^ defined as one of the 
columns A* with i G Np,, all columns with i G Np being multiples of each other. 
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(i ! ) Each component jl of C(A) is contained in a (spectral) class \i with p® = 1, and 
the support of each x^ of (i) consists of all indices in all classes that have access 
to /i. 

(ii) V©(A, 1) is generated by of (i), for jl ranging over all components of C(A). 

(iii) x^ of (i) are extremals in V (S (A,1). (Moreover, a?w are strongly linearly inde- 
pendent in the sense of \J5\.) 

To verify (i'), not explicitly stated in the references, use (i) and the path interpretation 
of A*. 

Vectors x^ of Theorem 12.111 are also called the fundamental eigenvectors of A, in max 
algebra. Applying a strict visualization scaling allows to get further details on these 
fundamental eigenvectors. 

Proposition 2.12 ([41J Proposition 4.1). Let A G M+ Xn have p®{A) = 1 and be strictly 
visualized. Then 

(i) For each component fx ofC(A), x^"' of Theorem \2.11\ can be canonically chosen as 
A* for any i G Njx, all columns with i G Njx being equal to each other. 

(ii) Xj® < 1 for alii. Moreover, supp(x^W) = iV^. 

Further we also deal with the (Minkowski) sum of the eigencones associated with all 
p G A(A). Defining it simultaneously in both algebras, 

(7) V*(A):= £ V(A,p) 

peA(A) 

The cone V T '{A) consists of all sums of vectors in different V(A,p). 

2.6. Notation table. For the reader's convenience we give the following notations table 
showing how various objects are denoted in nonnegative algebra, max algebra and when 
both algebras are considered simultaneously. 
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span e (A) 
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Perron root 


pi 


Pv 

• I' 


Pu 


Spectrum (excl. 0) 


A+(A) 




A(A) 


Eigencone 


v + (Ap + ) 




V(A,p) 


Sum of eigencones 


V*{A) 


V^{A) 


V*(A) 


Core 


core+(A) 


core® (A) 


core (A) 



Classes /i, Critical comp. jl, Classes /i, 

In the case of max algebra, we also have the critical graph C(A) (with related concepts 
and notation), and booleanizations A^\ x 1 ^ of matrices and vectors A,x with all entries 
not exceeding 1 (not used in nonnegative algebra). 

3. Periodic sequences of eigencones 

3.1. Spectral classes and critical components of matrix powers. The aim of this 
subsection is to demonstrate that access relations and spectral classes of all matrix powers 
are essentially the same, and that the case of an arbitrary eigenvalue reduces to the case 
of the principal eigenvalue. See Theorems 13.51 and 13.61 for main formulations. 

Results of this subsection can be derived, for the case of nonnegative algebra, from 
the classical Perron- Frobenius theory. Some of them can be traced back to the classical 
work of Frobenius |17j. To show this we now quote (in translation) some remarks on 
the very first page of [17] concerning the powers of an irreducible nonnegative matrix 
A. Frobenius defines (what we could call) the cyclicity or index of imprimitivity k of 
an irreducible S as the number of eigenvalues that lie on the spectral circle. He then 
remarks "If A is primitive, then every power of A is again primitive and a certain power 
and all subsequent powers are positive" . This is followed by "If A is imprimitive, then A m 
consists of d irreducible parts where d is the greatest common divisor of m and k. Further, 
A m is completely reducible. The characteristic functions of the components differ only in 
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the powers of the variable" (which provides a converse to the preceding assertion). And 
then "The matrix A k is the lowest power of A whose components are all primitive" . The 
three quotations cover Lemma [3.11 below in the case of nonnegative algebra. 

The reader is also referred to the monographs of Mine [28j , Berman-Plemmons [3] , Brua- 
ldi-Ryser [4] , and we will often cite the work of Tarn-Schneider |43] Section 4 containing 
all of our results in this section, in nonnegative algebra. 

Lemma 3.1 (cf. [3] Ch. 5, Ex. 6.9; [13] Lemma 4.5). Let A be irreducible with the unique 
eigenvalue p, let Q(A) have cyclicity o and k be a positive integer. 

(i) A k is a direct sum of gcd(k,a) irreducible blocks with eigenvalues p k . The set of 
all indices in these blocks is {1, . . . ,n} , and A k does not have eigenvalues other 
than p k . 

(ii) If k is a multiple of a, then the sets of indices in these blocks coincide with the 
cyclic classes ofQ(A). 

(hi) If supp(x) is a cyclic class ofQ(A), then supp(v4a;) is the previous cyclic class. 

Proof, (i): Assuming w.l.o.g. p = 1, let X = diag(a;) for a positive eigenvector x £ V(A, p) 
and consider B := X~ X AX which is stochastic (nonnegative algebra), or max-stochasticL 1 
(max algebra). All gcd(/c, a) irreducible isolated blocks of B k (Theorem l2.2l) are stochastic 
(or max-stochastic), hence they all have an eigenvector (1, . . . , 1) associated with the 
eigenvalue 1. If x £ V(A k ,p k ), then its subvectors corresponding to the irreducible blocks 
of A k are also eigenvectors of those blocks. Hence p = 1 is the only eigenvalue. 

(ii) : By Theorem 12.21 Q(A) splits into gcd(/c, a) = o components, and each of them 
contains exactly one cyclic class of Q{A). 

(iii) : Use the definition of cyclic classes and that each node has an ingoing edge. 

□ 

In both algebras, each class p of A with cyclicity a corresponds to an irreducible 
submatrix A^. It is easy to see that (A k ) flfl = (A^) k . Applying Lemma EH] to A w we 
see that p gives rise to gcd(/c, a) classes in A k , which are said to be derived from their 



1 such that bij = 1 holds for all 



TWO CORES OF A NONNEGATIVE MATRIX 17 

common ancestor ii. The classes of A k and A 1 derived from the common ancestor will be 
called related. Note that this is an equivalence relation on the set of classes of all powers 
of A. 

In max algebra, with A assumed to be strictly visualized, each component jx of C(A) 
with cyclicity a corresponds to an irreducible submatrix (recall that A^ := A N . Njl ). 
It is easy to see that (A® k )f~ = (A l ~l)® k . Applying Lemma Oto Af- we see that \x gives 
rise to gcd(A;, a) critical components in A® k , which are also said to be derived from their 
common ancestor ft. 

Evidently, a class of A h , resp. a component of C(A® k ), is derived from a class of A or a 
component of C(A), if and only if its index set is contained in the index set of the latter 
class (or the component) 

Lemma 3.2. Both in max algebra and in nonnegative linear algebra, sets of indices in 
the classes of A k with nonzero Perron roots coincide for all k. 

Proof. In both algebras, an index belongs to a class with nonzero Perron root if and only 
if the associated graph contains a cycle with a nonzero weight traversing the node with 
that index. This property is invariant under taking matrix powers, hence the claim. □ 

Theorem 3.3 (cf. [6] Theorem 8.2.6, [Sj Theorem 2.3). Let A G M+ Xn . 

(i) Sets of critical indices of A® k coincide for all k. 

(ii) Each component fi of C(A) with cyclicity o is associated with a set of gcd(k, a) 
components of C(A® k ) derived from it, and each component of C(A® k ) is derived 
from a component of C(A). If k is a multiple of a, then the index sets in the 
derived components are the cyclic classes of jl. 

(iii) If A is strictly visualized, Xj < 1 and supp(xW) is a cyclic class of fx, then 
supp((Ar)W) is the previous cyclic class of fx. 

Proof, (i): Use that C{A® k ) = (C(A)) k . (ii): Apply Theorem IO to a component p, of 
C(A). (iii): Use the definition of cyclic classes and that each node in jl has an ingoing 
edge. □ 
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We observe that access relations in matrix powers are similar. This has identical proof 
in max algebra and nonnegative algebra. 

Lemma 3.4. For all k, I > 1 and p > 0, an index i G {l,...,n} accesses (resp. is 
accessed by) a class with Perron root p k in A k if and only if it accesses (resp. is accessed 
by) a related class with Perron root p l in A 1 . 

Proof. We treat only the situation where i accesses the class, the remaining case is anal- 
ogous. Our argument goes from a class of A k back to the ancestor class of A, and then 
forward to a related class of A 1 . 

We deduce from Lemma 13.11 and Lemma 13.21 that the index set of each class of A k 
with Perron root p k is contained in the ancestor class of A with Perron root p. Then, i 
accessing a block in A k implies i accessing its ancestor in A. Since p > 0, this ancestor 
class is nontrivial, so the access path can be extended to have a length divisible by I, by 
means of a path contained in the ancestor class. By Lemma I3TT] the ancestor decomposes 
in A 1 into several classes with the common Perron root p l , and i accesses one of them. □ 

Both in max algebra and in the nonnegative linear algebra, the spectrum of a nonneg- 
ative matrix is determined by 1) Perron roots of the blocks, 2) access relations to the 
blocks with a particular Perron root. 

Theorem 3.5 ([43J, Corollary 4.6). Let A G M" xn . 

(i) For any p G R + , the sets of all indices belonging to the spectral classes of A k 
associated with p k are the same for all k. 

(ii) A(A k ) = {p k ; p G A(A)}. 

(iii) For each spectral class p of A with cyclicity a there are gcdik, a) spectral classes of 
A k derived from it. Conversely, each spectral class of A k is derived from a spectral 
class of A. If k is a multiple of a then the index sets of derived spectral classes 
are the cyclic classes of p. 

Proof, (i): Let i G {1, . . . , n}, and denote by p(A, i) the Perron root of the class of A to 
which % belongs. Lemmas 13.11 and 13.21 imply that p(A k ,i) = p k (A,i) for all k. Hence for 
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any p G A(A) the sets of indices i with p(A k , i) = p k and p(A l , i) = p l are the same for all 
k, I > 1, and it remains to show that the same holds if we restrict to spectral classes. 

It suffices to consider the case when i belongs to a nontrivial class (in all powers). An 
index i does not belong to a spectral class of A k if it is accessed in A k by a class p for 
which p M = p k > p k in max algebra, resp. p M = p k > p k in the nonnegative algebra. Here 
p p is the Perron root of class p in A k , and p is the Perron root of the ancestor of p. By 
Lemma T3.4[ i is accessed in A 1 by a class v related to p, and hence such that p v = p l > p l 
in max algebra, resp. p v = p l > p l in the nonnegative algebra. Here p v is the Perron root 
of class v of A 1 , and p is still the Perron root of the common ancestor of p and v. Hence i 
does not belong to any spectral class of A 1 . Thus the set of indices in all spectral classes 
of A k with Perron root p k coincides with the set of indices in all spectral classes of A 1 
with eigenvalue p l for all k, I > 1. 

(ii) : We have p(A k ,i) = p k (A,i) for all i (see the beginning of the proof of (i)), and by 
(i) the sets of all indices in spectral classes are the same for all powers. 

(iii) : Lemma [3.11 shows that the set of all indices in the gcd(fc, a) classes of A k derived 
from some spectral class of A with cyclicity a is the same as the set of all indices in that 
spectral class. By (i) the sets of all indices in spectral classes coincide, in particular the 
derived classes are spectral, and there is no room for any spectral classes of A k other than 
derived from the spectral classes of A. □ 

Now we show that as in the case of eigencones of a matrix, working with V(A k ,p k ) we 
can assume that p = 1 is the principal eigenvalue of A, and hence of all A k . 

Theorem 3.6. Let k > 1 and p G A(A). 

(i) The set of all indices having access to the spectral classes of A k with the eigenvalue 
p k equals M p , for each k. 

(ii) (A k ) MpMp =p k (A p ) k MpMp . 

(iii) V(A k ,p k ) = V((A p ) k A). 
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Proof, (i): Apply Theorem 13.51 part (i) and Lemma 13.41 (ii): Use that M p is initial in 
g{A). (iii): By Proposition EH we have V{A k ,p k ) = V{A k pk , 1) where, instead of ©, 

A k k := p~ k ( ° ° | , and 

(8) V° a *m) 

M k := {i; i -> v, v is (A k , p fc )-spectral} 
By part (i) M k = M p , hence A k k = (A p ) k and the claim follows. □ 

3.2. Periodicity of eigencones. Here we investigate the periodicity of eigencones. Un- 
like in the previous subsection, the proof of periodicity will be different for the cases of 
max algebra and nonnegative algebra. The periods of eigencone sequences in max algebra 
and in nonnegative linear algebra are also in general different, for the same nonnegative 
matrix (see Sect. 4.4 for an example). 
Observe that in both algebras 

(9) k divides I V(A k , p k ) C V{A\ p l ) Vp e A(A) 

We start with the periodicity of eigencones in max algebra. 

Theorem 3.7. Let p G A©(A) and a be the cyclicity of the corresponding critical graph. 
Then V @ {A® l ,p l ) = U ffi (A®( /+(T ), p /+CT ) and V & (A & , p l ) C V^A® , p ff ) for all I. 



Proof. We can assume w.l.o.g. that p = 1, and using Theorem l3.6l that this is the greatest 
eigenvalue of A. 

By Theorem I2.5[ columns of A® r with indices in N C (A) are periodic at r > T C (A). 
Recall that T C (A) is not less than the periodicity threshold T(C(A)). By Theorem 12.11 
part (ii), (C(A)) ta consists of complete graphs for ta > T(C(A)), in particular, it contains 
loops for all i G N C (A). Hence 

af tu) = 1 ViG N C (A), ta > T C (A), 

and 

*T ta) > = < Vz G N C (A), Vfc, I, Via > T C (A), 
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or, in terms of columns of matrix powers, 

A a(!+t*) > A & Vi e Nc (A), V/, \/ta > T C {A) 

Multiplying this inequality repeatedly by A® 1 we obtain A ®.^ kl+ta ^ > j\®( kl ) f or all A; > 1, 
or Af k{l+ta)) > Af kl) for all k > 1. Hence we obtain 

(10) > (^4®')* Vi G iV c (A), V/, Via > T C (A) 

On the other hand, using the ultimate periodicity of critical columns we have 

= 0{Af ; s = fcZ(mod(7), fc > 1, s > T C (A)} 
for all / and all ter > T C (A), while generally 

(A®')* > 0{Af ; s = fcZ(modo-), k > 1, s > T C (A)} V/, 
implying the opposite of ffTUl) . It follows that 

(11) ( A m+t<r)}*. = Vi e iVc (^) ) v/, vta > t c (a), 

therefore (A®(' +<7 ))* = (A®^*^)* = (A®^ 1 ^ = (A^)* for all critical indices % and all 
I. Since , 1) is generated by the critical columns of (A® 1 )*, and the critical indices 

of A® 1 are N C (A) by TheoremESJ the periodicity V e (A® l ,p l ) = V^A®^, p l+a ) follows. 
Using this and © we obtain V @ (A & , p l ) C V®(A®^, p i<T ) = V r e (A® a , p CT ) for each Z and 

p e A©(A). □ 

We now show a nonnegative version of Theorem 13.71 

Theorem 3.8. Let p G A + (v4) and cr fre i/ie 1cm of the cyclicities of all the corresponding 
spectral blocks. Then V + {A xl ,p l ) = V + {A x( - l+C7 \ p l+a ) and V + {A xl ,p l ) C V + (A Xa ,p a ) for 
all I. 

Proof. Assume w.l.o.g. that p = 1, and using Theorem 13. 6[ that this is the greatest 
eigenvalue of A and all final classes (and hence only them) have Perron root p. Final 
classes of A xl are derived from the final classes of A; they (and no other classes) have 
Perron root p l . By Theorem 13.51 and Corollary 12. 3l the spectral classes of A xl and A x ( l+ta \ 
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for any t > 0, have the same index sets, which we denote by Ni, . . . , N m (assuming that 
their number is m > 1). 

By the Frobenius- Victory Theorem 12.101 the cone V + (A xl , 1) is generated by m eigen- 
vectors x^\ . . . with a support condition, from which we infer that the subvectors 
Xpfi (i.e., Xj!f) are positive, while xi^ (i.e., xj^J) are zero for all p, ^ v from 1 to m, since 
the different spectral classes by (J5J) do not access each other, in the nonnegative linear 
algebra. Analogously the cone V + {A x ^ l+t<T \ 1) is generated by m eigenvectors . . . , y( m > 
such that the subvectors are positive, while yi = for all p ^ v from 1 to m, 

Let I = a. As V + {A Xa , 1) C V + (A^ ta \ 1), each x^ is a nonnegative linear combination 
of y^\ . . . , y( m \ and this implies x^> = y^ for all p = 1, . . . , m. Hence V + (A x ^ ta \ 1) = 
V+(A X<7 ,1) for all t > 0. 

We also obtain V + {A X \ 1) C y + (A x ^, 1) = 1/ + (A XCT ,1) for all Z. Thus V + {A X \ 1) C 
^+ (A X (* CT ), 1), and therefore V+(A xi , 1) C V + (A x( - l+t ^ , 1). Now if V + (A xi , 1), resp. 
V+(A x ^ + '°"), 1) are generated by r 1 ',...,^™', resp. y^\ . . . , y^ described above and 
each x^ is a nonnegative linear combination of . . . , f/ m \ this again implies x^ = y^ 
for all p = 1, . . . , m, and V + (A x( - l+ta \ 1) = V + (A xl , 1) for all t > and all I. The inclusion 
V + (A xl , 1) C V+(A Xa , 1) was shown above. 

□ 

Based on Theorems 13.71 and I3.8[ let us introduce the following notations, which can be 
seen as definition of cyclicity in nonnegative algebra and in max algebra. 

1. Let <j p be the the 1cm of all cyclicities of spectral classes with p G A+(A) (non- 
negative algebra), or the 1cm of all cyclicities of critical graphs associated with 
p G A ffl (A) (max algebra). 

2. Let o"a be the 1cm of all a p defined in 1. 
Note the following version of (Q 

(12) k divides Z => V*(A k ) C V E {A l ) 

for the sum of all eigencones V T '(A k ) = XlpeA(.4) V(A k ,p k ). Summarizing Theorems 13.71 
and 13.81 and using a p and a\ we obtain the following formulation (now uniting both 
algebras) . 
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Theorem 3.9. Let A G M+ Xn . For all I > 1 and p G A(A), 

(i) V(A z ,/0 Z ) = and F(A',p') C 1/(A^,p CT "), 

(ii) V^{A l ) = 1/ E (A' +CTA ) and V^{A l ) C \/ s (A CTA ). 

3.3. Inclusion and divisibility. We now show that the inclusion relations between the 
eigencones of different powers of a matrix, in both algebras, strictly follow divisibility of 
exponents of matrix powers with respect to o p and a\. 
Here is another consequence of Theorems 13.71 and 13.81 

Lemma 3.10. Let k,l>l and p G A(A). 

(i) V(A k ,p k ) = V(A 9cd ^ a P' k \ p 9cd{o P ,k)^ and v^(A k ) = v Y; (A 9cd ^ A ' k) ). 

(ii) gcd(k,a p ) = gcd(l,a p ) implies V(A k ,p k ) = V(A l ,p l ), and gcd(k,a\) = gcdi^l.a^) 
implies V E (A k ) = V^(A l ). 

Proof, (i): Let a := o~ p , and s :=gcd(fc,cr). If s = a then k is a multiple of a and 
V(A k ,p k ) = V(A s ,p s ) by Theorems 13.71 and 13.81 there is nothing to prove. Otherwise, 
since s divides k, we have V(A S , p s ) C V(A k , p k ). In view of the periodicity (Theorems 13.71 
andEE}, it suffices to find t such that V(A k ,p k ) C V(A s+ta , p s+ta ). For this, observe that 
s + to is a multiple of s = gcd(fc,<r). By Lemma 12.71 (the coin problem of Frobenius), 
for big enough t it can be expressed as t\k + t 2 cr where t±,t2 > 0. Moreover t± ^ 0, for 
otherwise we have s = o. Then we obtain 

V(A k ,p k ) C V(A hk ,p hk ) = v(A hk+t2a ,p hk+t2a ) 

= V{A s+ta ,p s+ta ) = V{A s ,p s ), 

and the first part of the claim follows. The second part is obtained similarly, using 
Theorem 13.91 part (ii) instead of Theorems 13.71 and 13.81 

(ii) follows from (i). □ 

Proposition 3.11. For each k > 1, the set of extremals ofV T '(A k ) is the union of the 
sets of extremals ofV(A k ,p k ) for p G A(A). 



Proof. We can assume w.l.o.g. that k = 1, due to Theorem 13.51 (ii). 
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1. As V T '{A) is the sum of V(A,p) for p G A(A), it is generated by the extremals of 
V(A, p) for p G A (A). Hence each extremal of V^ E (A) is an extremal of V(A, p) for some 
p G A(v4). 

2. Let x G V(-A, p^), for some spectral class p, be extremal. Assume w.l.o.g. that 
Pn — 1, and by contradiction that there exist vectors y K , all of them extremal eigenvectors, 
such that x = Y2 K V* If there is y K associated with an eigenvalue p v > 1, then applying A* 
we obtain x = (puYy K + ■ ■ -, which is impossible at large enough t. So p v < 1. With this 
in mind, if there is y K associated with p v < 1, then 1) in nonnegative algebra we obtain 
Ac > A y K , a contradiction; 2) in max algebra, all nonzero entries of Ay K go below 
the corresponding entries of x meaning that y K is redundant. Thus we are left only with 
y K associated with p v = 1, which is a contradiction: an extremal x G V(A, 1) appears as 
a sum of other extremals of V(A, 1) not proportional to x. □ 

Theorem 3.12. Let A G IR™ xn and a be either the cyclicities of spectral classes of A (non- 
negative algebra) or the cyclicities of critical components of A (max algebra) . The 

following are equivalent for all positive k, I: 

(i) gcd(k, a) divides gcd(l, a) for all cyclicities a; 

(ii) gcd(k,a p ) divides gcd(l,a p ) for all p G A(A); 

(iii) gcd(k,o-^) divides gcd{l,ax); 

(iv) V(A k ,p k ) C V(A\p l ) for all p G A(A) and 

(v) V^(A k ) C V^(A l ). 

Proof. (i)=^(ii)=^(iii) follow from elementary number theory. (ii)=>- (iv) follows from 
and Lemma [3]T0] part (i), (iii)=>(v) follows from Lemma fT2l and Lemma [3 . 1 01 part (i). Note 
that (ii)=> (iv) and (iii)=^(v) are essentially based on Theorems 13.71 and 13.81 (iv)=^(v) is 
trivial. It only remains to show that (v)=> (i). 

(v)=> (i): In both algebras, take an extremal x G V(A k ,p k ). As V T '(A k ) C \/ s (A), 
this vector can be represented as x = where are extremals of V rS (v4'). There 

is a unique spectral class p of A to which all indices in supp(x) have access, and since 
supp y l C supp x, we are restricted to the submatrix Ajj where J is the set of all indices 
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accessing p in A. In other words, we can assume w.l.o.g. that p is the only final class in 
A, hence p is the greatest eigenvalue, and p = 1. Note that suppx PI ^ 0. 
In nonnegative algebra, restricting the equality x = Yli V % to we obtain 

(13) supp(a; M ) = (Jsupp(j£). 

i 

If supp(?/^) is non-empty, then y % is associated with a spectral class of A xl whose nodes 
are in iV M . Theorem 12.101 implies that supp(y^) consists of all indices in a class of A* 1 . 
Taking as x any extremal eigenvector of A xk with supp x n ^ 0, fflBl shows that each 
class of A** (corresponding to x) splits into several classes of A* (corresponding to y l ). 
By Corollary 12.31 this is only possible when gcd ( k, a) divides gcd(/,cr), where a is the 
cyclicity of the spectral class p. 

In max algebra, since p = 1, we can w.l.o.g. assume that A is strictly visualized. In 
this case A and x have all coordinates not exceeding 1. Recall that x^ is the Boolean 
vector defined by xf 1 = 1 <^ Xi = 1. Vector x corresponds to a unique critical component 
p of C(A) with the node set Njx- Then instead of ( JTBl) we obtain 

(14) rrW=0y^ =► supp^ 1 ) = [j supp^ 11 ), 

i i 

where supp (a; W) = supp(x^) by Proposition 12.121 and Theorem 13.31 and hence also 
supp (y l M) = supp(y~ ). If supp(y~ ) is non-empty then also supp(?/]y ) is non-empty 
so that y % is associated with the eigenvalue 1. As y % is extremal, Proposition 12.121 implies 
that supp(y~^) consists of all indices in a class of (A [ }1)® 1 . Taking as x any extremal eigen- 
vector of A® k with supp(xW) Pi Njx ^ 0, ffT4l) shows that each class of (^4- i 1 i)® fc splits into 
several classes of (A^~~~)® 1 . By Corollary 12.31 this is only possible when gcd (A;, a) divides 
gcd(/, a), where a is the cyclicity of the critical component p. □ 

Let us also formulate the following version restricted to some p G A(A). 

Theorem 3.13. Let A G M™ xn ; and let a be either the cyclicities of spectral classes (non- 
negative algebra,) or the cyclicities of critical components (max algebra^ associated 
with some p G A(A). The following are equivalent for all positive k, I: 

(i) gcd(k, a) divides gcd(l, a) for all a; 
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(ii) gcd(k,a p ) divides gcd(l, a p); 

(iii) V(A k ,p k ) C V(A l ,p l ). 

4. Core and eigenvectors 

4.1. Basics. In this section we investigate the core of a nonnegative matrix defined as 
the intersection of the cones generated by its powers: 

(15) core(A) := span(A i ) 

In the main argument, we consider the cases of max algebra and nonnegative algebra 
simultaneously. 

One of the most elementary and useful properties of the intersection on the r.h.s. of ( fl5|) 
is that, actually, 

(16) span(A) D span(A 2 ) D span(A 3 ) D ... 
Slightly generalizing an argument of Pullman [33] we will show that 

(17) core(A) = ^V s (A fc ) = ^ V(A k ,p k ) 

k>l fc>l,peA(A) 

also in max algebra. 

Note that the following inclusion is almost immediate. 

Lemma 4.1. Efc>i^ S (^ fc ) Q core(A). 

Proof, x G V(A fc ,p) implies that A k x = px and hence x = p~ l A kt x for all t > 1 (using 
the invertibility of multiplication). Hence x G f\>i spanv4 fc * = f\>i span(74*). □ 

So it remains to show the opposite inclusion 

(18) core(A) cY^V s (A k ). 

k>\ 

Let us treat the trivial case p(A) = 0, i.e., A(A) = 0. There are only trivial classes in 
the Frobenius normal form, and G(A) is acyclic. This implies A k = for some k > 1. In 
this case core(A) = {0}, all V(A k ,p k ) with p > are {0}, so (EE7j) is the trivial "draw" 
{0} = {0}. 



TWO CORES OF A NONNEGATIVE MATRIX 27 

4.2. Max algebra: cases of ultimate periodicity. In max algebra, unlike the non- 
negative algebra, there are wide classes of matrices for which (ITS]) and (ITTj) follow almost 
immediately. We list some of them below. 

51 : Irreducible matrices. 

52 ■ Ultimately periodic matrices. This is when we have A® t+a = p a A® t for some a and p, 
and all sufficiently large t. As shown by Molnarova-Pribis [22], this happens if and only 
if the Perron roots of all nontrivial classes of A equal = p. 

53 : Robust matrices. For any vector x G IR+ the orbit {A® 1 £g> 2} hits an eigenvector 
of A, implying that the whole remaining part of the orbit consists of multiples of that 
eigenvector. The notion of robustness was introduced and studied in [9]. 

154 : Orbit periodic matrices: For any vector x G IR+ the orbit {A® l ®x} hits an eigenvector 
of A® a<c for some a x , implying that the remaining part of the orbit is periodic (with some 
growth rate). See [ID] Section 7 for characterization. 

155 : Column periodic matrices. This is when for all i we have (A lX, ^ +f7 ^).j = Pi'A®* for all 
large enough t and some pi and crj. 

Observe that Si C S 2 C S4 C S 5 and <S 3 C £4. Indeed, Si C (S 2 is the Cyclicity 
Theorem 12.41 For the inclusion S2 Q S4 observe that, if A is ultimately periodic then 

A ®(t+a) = p a A ®t and hence A ®(t+a) ^ % = p a A ®t g, ^ holdg for gjj ^ ^ &nd &U big 

enough t. Observe that S3 is a special case of £4, which is a special case of 1S5 since the 
columns of matrix powers can be considered as orbits of the unit vectors. 

To see that (|T8|) holds in all these cases, note that in the column periodic case all 
column sequences {A 1 A- end up with periodically repeating eigenvectors of A® ai , which 
implies that span^A®*) C fc>1 V^(A® k ) C core®(A) and hence span^A®*) = core ffi (A) 
for all large enough t. Thus, finite stabilization of the core occurs in all these classes. A 
necessary and sufficient condition for this finite stabilization is described in |10j . 

4.3. Core: a general argument. The original argument of Pullman [33] Section 2 used 
the separation of a point from a closed convex cone by an open homogeneous halfspace 
(that contains the cone and does not contain the point). 
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In the case of max algebra, Nitica and Singer [3T1 [32] showed that at each point x G 
there are at most n maximal max-cones not containing this point. These conic semispaces, 
used to separate x from any max cone not containing x, turn out to be open. Hence they 
can be used in the max version of Pullman's argument. 

However, for the sake of a simultaneous proof we will exploit the following analytic 
argument instead of separation. By B(x, e) we denote the intersection of the intersection 
of the open ball centered at x G with radius e with M™. In what follows up to 
Theorem 14.61 we consider both algebras simultaneously. 

Lemma 4.2. Let x 1 , . . . , x m G R" be nonzero and let z ^ span(x 1 , . . . , x m ). Then there 
exists e > such that z span(i?(x 1 , e) U . . . U B(x m , e)). 

Proof. By contradiction assume that for each e there exist points y l {e) G B(x\e) and 
scalars /ij(e) such that 

m 

(19) z = J>(ey(e). 

i=l 

Since y l (e) — > x l as e — > and x l are nonzero, we can assume that y l (e) are bounded 
from below by nonzero vectors v l , and then z > Y^Li l Ji i( e ) vl f° r a ^ e > implying that 
/ii(e) are uniformly bounded from above. By compactness we can assume that Hi(e) 
converge to some /ij G M+, and then (fT9!) implies by continuity that z = X^i/ 1 * 2 ' 4 ? a 
contradiction. □ 

Theorem 4.3 Q33J Theorem 2.1). Assume that {K{\ for I > 1, is a sequence of cones in 
1R™ such that Ki + \ C K\ for all I, and each of them generated by no more than k nonzero 
vectors. Then the intersection K = D^Ki is also generated by no more than k vectors. 

Proof. Let K\ = span(y n , . . . , y lk ), and consider the sequences of normalized vectors 
{y h /\ \y h \ I? ; ^ > 1} for i = where := maxiij (or any other norm). As 

the set {u; \\u\\ = 1} is compact, we can find a subsequence l t such that y'*V||y'* l ||, 
as t — > oo, converge to a finite vector u 1 , which is nonzero since = 1. We claim 
that K = span(w 1 , . . . ,u k ). Assume to the contrary that there is z G K that is not in 
span(w 1 , . . . ,u k ). Then by Lemma 14.21 there exists e > such that z ^ span(S(M 1 , e) U 
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. . . U B(u k ,e)). Since y ln converge to it* as t — > oo, we have y ln G B(u l ,e) for t large 
enough, and 

span(i/ u , . . . , y ltk ) C span^n 1 , e) U . . . U 5(« fe , e)) 

But ^ belongs to Ki t = span(y Ztl , . . . , i/' ) since it belongs to the intersection of all these 
cones, a contradiction. □ 

Theorem 14.31 applies to the sequence {span(y4*), t > 1}, so core (.A) is generated by no 
more than n vectors. 

Proposition 4.4 ([33], Lemma 2.3). The mapping induced by A on its core is a surjection. 

Proof. First note that A does induce a mapping on its core. If z G core (A) then for each t 
there exists x* such that z = A t x t . Hence Az = A t+1 x l , so Az G C\ t >2 span A 1 = core(A). 

Next, let m be such that A m has the greatest number of zero columns (we assume that 
A is not nilpotent; recall that a zero column in A k remains zero in all subsequent powers). 
If z = A t x t for t > m + 1, we also can represent it as y4 m+1 it*, where u f := A t ~ m ~ 1 x t . The 
indices of it* corresponding to the nonzero columns of A m+1 are bounded since A m+1 u t = z. 
So we can assume that the sequence of subvectors of it* with these components converges. 
Then the sequence y* := A m u l also converges, since the indices of nonzero columns of A m 
coincide with those of A m+l , which are the indices of the converging subvectors of it*. Let 
y be the limit of y l . Since y s = A s ~ l x s are in span(A*) for all s > t, and since span(A*) are 
closed, we obtain y G span(v4*) for all t. Thus we found y G core(A) satisfying Ay = z. □ 

Theorem 14.31 and Proposition 14.41 show that the core is generated by finitely many 
vectors in and that the mapping induced by A on its core is "onto" . 

Now we use that a finitely generated cone in the nonnegative orthant (and more gen- 
erally, closed cone) is generated by its extremals both in nonnegative algebra and in max 
algebra, see ff\ |4"4"]. 

Proposition 4.5 ([33], Theorem 2.2). The mapping induced by A on the extremal gen- 
erators of its core is a permutation (i.e., a bijection). 

Proof. Let core(v4) = span(it 1 , . . . , it fe ) where it 1 , . . . ,u k are extremals of the core. Suppose 
that x J is a preimage of it 7 in the core, that is, Ax^ = it 7 for some x J G coie(A), j = 
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1, . . . , k. Then x J = £)i=i a i u% f° r some nonnegative coefficients a%, . . . , a^, and m j = 
5^*=i O-iAu 1 . Since u J is extremal, it follows that u J is proportional to Au % for some i. 
Thus we have a well-defined inverse mapping on the set of extremals, which has to be 
injective. Since the number of extremal rays is finite, this inverse mapping is also bijective 
on them. □ 



We are now ready to show that (ITTI) holds in both algebras. Below we use the notation 
o"a and later a p (see before Theorem 13.91) 

Theorem 4.6. Let A G M+ Xn . Then 

core(A) = V ( Ak 'P k ) = V s {A aA ). 

fc>l,p€A(A) 

Proof. Proposition 14. 51 implies that all extremals of core(^4) are eigenvectors of some power 
of A, for which we take the 1cm of the lengths of all cycles of the permutation induced 
by A on the extremals. Hence core (A) is a subcone of the sum of all eigencones of all 
powers of A, which is the inclusion of f ll8p . Combining this with the converse inclusion 
of Lemma [4.11 we obtain that core(A) is precisely the sum of all eigencones of all powers 
of A ( TT~T|) . The last part of the equality now follows from the periodicity of eigencones 
described in Theorem 13.91 □ 

The nonnegative part of Theorem 14.61 can be found in Tarn-Schneider ( |43] . Th. 4.2 
part (hi)), and we now prove an extension of ([HI], Th. 4.7) covering both algebras. A 
vector x G is called normalized if maxxj = 1. 

Theorem 4.7. Let A G M+ Xn . 

(i) The set of extremals of core(A) is the union of the sets of extremals of V(A a , p u ) 
for all p G A(A), with a = a p . 

(ii) In nonnegative algebra, each spectral class p, with cyclicity a p corresponds to a 
set of distinct a p normalized extremals of core + (A) , with all indices in the support 
having access to p. In max algebra, each critical component p with cyclicity 
a p associated with some p G A^(A) corresponds to a set of distinct normalized 
extremals x of core® (A), which are (normalized) columns of (Ap p )* with indices in 
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(iii) Each set of extremals described in (ii) forms a simple cycle under the action of A. 

(iv) There are no normalized extremals other than those described in (ii). In non- 
negative algebra, the total number of normalized extremals equals the sum of 
cyclicities of all spectral classes of A. In max algebra, the total number of 
normalized extremals equals the sum of cyclicities of all critical components of A. 

Proof, (i) follows from Theorem 14.61 and Proposition 13.111 

For the proof of (ii) and (iii) we can fix p = p p G A (A), assume A = A p (using 
Theorem I3.6p and a := a p (defined differently in two algebras, see Theorem 13 .9p . In max 
algebra, we also assume that A is strictly visualized. 

(ii) In nonnegative algebra, observe that by Theorem I3.5[ each spectral class p of A 
gives rise to a p spectral classes in A Xa (note that a p divides a). According to Frobenius- 
Victory Theorem 12. 101 and Lemma \3. 41 these classes, whose node sets are cyclic classes of 
p, give rise to normalized extremals of V + (A XeT , 1) with the given support condition. 

(iii) : Let x be an extremal described above. Then supp(x)nA^ is a cyclic class of p and 
supp(Ar) D Nfj, is the previous cyclic class of p, by Lemma I3TT1 part (iii). It can be checked 
that all indices in supp^rc) also have access to p. By Proposition 14. 5[ Ax is an extremal 
of core + (A), and hence an extremal of V + (A Xcr , 1). Theorem 12. 101 and Lemma \3 . 1 1 identify 
Ax with the extremal associated with the previous cyclic class of p. 

Vectors x, Ax, . . . , A Xa ^~ 1 x are distinct since the intersections of their supports with 
N p are disjoint, so they are exactly the set of extremals associated with p. Note that 
A Xa > J -x = x, as supp(A Xcrfl a;) n iV M = supp(x) D N p , and both vectors are extremals of 
V+{A X °A). 

(ii) : In max algebra, observe that by Theorem 13.31 each component p of C(A) gives 
rise to a p components C(A® a ), whose node sets are the cyclic classes of p (note that 
o-jx divides a). These components correspond to a p columns of (A 00 ")* with indices in 
different cyclic classes of p, which are by Theorem 13.31 the node sets of components of 
C(A® a ). By Theorem 12.111 and Proposition 12. 12[ these columns of (A 8 " 7 )* are normalized 
extremals of V r e (A® a , 1). 

(iii) : Let x be an extremal described above. By Proposition 12.121 and Theorem 13.31 
supp(xW) is a cyclic class of p, and supp((A <g) x)W) is the previous cyclic class of p. By 
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Proposition 14. 5\ A £g> x is an extremal of core e (A), and hence an extremal of V§(A® a , 1). 
Proposition 12.121 and Theorem 13.31 identify A <g> x with the extremal associated with the 
previous cyclic class of fx. 



Vectors x, A 



x. 



A®<Tji i x are distinct since their booleanizations 



x [l \ (A <g> x) 



[i] 



are distinct, so they are exactly the set of extremals associated with fx. Note that A® a ^x = 
x, as (A^^^x)^ = a;M and both vectors are extremals of V <S) (A®' 7 , 1). 

(iv): In both algebras, Theorem 13.51 shows that there are no spectral classes of A a 
other than the ones derived from the spectral classes of A. In nonnegative algebra, 
this shows that there are no extremals other than described in (ii). In max algebra, 
on top of that, Theorem 13.31 shows that there are no components of C(A® a ) other than 
the ones derived from the components C(A p ), for each p £ A 9 (A), hence there are no ex- 
tremals other than described in (ii). In both algebras, it remains to count the extremals 
described in (ii). □ 

4.4. Examples. We consider two examples of reducible nonnegative matrices, examining 
their core in max algebra and in nonnegative linear algebra. 
Example 1. Take 



(20) 



A 



/0.1206 











\ 


0.5895 


0.2904 


1 


0.8797 


0.4253 


0.2262 


0.6171 


0.3439 


1 


0.3127 


0.3846 


0.2653 


0.5841 


0.2607 


1 


\0.5830 


1 


0.1078 


0.5944 


0.1788/ 



A has two classes with node sets {1} and {2, 3, 4, 5}. Both in max algebra and in nonnega- 
tive linear algebra, the only spectral class arises from M := {2, 3, 4, 5}. The max-algebraic 
Perron root of this class is p®{A) = 1, and the critical graph consists of just one cycle 
2->3->4->5->l. 
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The eigencones Vq(A, 1), V^A® 2 , 1), V^(A® 3 , 1) and K ffi (A® 4 , 1) are generated by the 
last four columns of the Kleene stars A*, (A® 2 )% (A® 3 )*, (A® 4 )*. Namely, 

Ve(A, 1) = V®(A® 3 , 1) = span ffi {(0 1111)}, 
V e (A m , 1) = span e {(0, 1, 0.8797, 1, 0.8797), (0, 0.8797, 1, 0.8797, 1)}, 
V e (A m , 1) = span e {(0, 1, 0.6807, 0.7738, 0.8797), (0, 0.8797, 1, 0.6807, 0.7738), 
(0, 0.7738, 0.8797, 1, 0.6807), (0, 0.6807, 0.7738, 0.8797, 1)} 

By Theorem I4.6[ core^(A) is equal to V^A® 4 , 1). Computing the max-algebraic powers 
of A we see that the sequence of submatrices Afj M becomes periodic after t = 10. In 
particular, 

fa \ 

0.4511 0.7738 0.6807 1 0.8797 
(21) A® 10 = 0.5128 0.8797 0.7738 0.6807 1 

0.5830 1 0.8797 0.7738 0.6807 

\0.5895 0.6807 1 0.8797 0.7738/ 

where < a < 0.0001. Observe that the last four columns are precisely the ones that 
generate V^(A® 4 , 1). Moreover, if a was then the first column would be the following 
max-combination of the last four columns: 



-"--2 ^D* u 51 -3 ^ u 21 -4 ^ 31 -5 



On one hand, the first column cannot be a max-combination of the last four columns 
since a > 0. On the other hand, a — > as t — > oo ensuring that the first column belongs 
to the core "in the limit" . 

Figure [1] gives a symbolic illustration of what is going on in this example. 

In nonnegative algebra, the block Amm with M = {2, 3, 4, 5} is also the only spectral 
block . Its Perron root is approximately p + {A) = 2.2101, and the corresponding eigencone 
is 



V+(A,p + (A)) = span + {(0, 0.5750, 0.5107, 0.4593, 0.4445)}. 
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Figure 1. The spans of matrix powers (upper curve) and the periodic 
sequence of their eigencones (lower graph) in Example 1 (max algebra). 

Taking the usual powers of (A/p + (A)) we see that 

fa \ 

0.2457 0.2752 0.2711 0.3453 0.2693 

0.2182 0.2444 0.2408 0.3067 0.2392 

0.1963 0.2198 0.2165 0.2759 0.2151 

\0.1899 0.2127 0.2096 0.2670 0.2082/ 



(A/P + (A)) 



xl2 



where < a < 0.0001, and that the first four digits of all entries in all higher powers 
are the same. It can be verified that the submatrix (A / p + (A))^1 { is, approximately, the 
outer product of the Perron eigenvector with itself, while the first column is also almost 
proportional to it. 
Example 2. Take 



(22) 



A 



( 1 \ 

10 

0.6718 0.2240 0.5805 0.1868 

\0.6951 0.6678 0.4753 0.3735/ 
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This matrix has two classes p, and v with index sets {1, 2} and {3, 4}, and both classes 
are spectral, in both algebras. In max algebra p® = 1 and p® = a 33 < 1. The eigencones 
of matrix powers associated with p® — 1 are 

V®{A, 1) = span e {(l, 1, 0.6718, 0.6951)}, 
V & (A m , 1) = span ffi {(l, 0, 0.3900, 0.6678), (0, 1, 0.6718, 0.6951)}, 

and the eigencone associated with p® is generated by the third column of the matrix: 

V@{A,p®) = span e {(0, 0, 0.5805, 0.4753)}. 

By Theorem 14.6} core®(A) is equal to the (max-algebraic) sum of V^(A® 2 , 1) and 
Vq(A, p®). To this end, observe that already in the second max-algebraic power 

/ 1 \ 

10 

0.3900 0.6718 0.3370 0.1084 

\0.6678 0.6951 0.2759 0.1395/ 



(23) 



A 



®2 



the first two columns are the generators of V^(A® 2 , 1). However, the last column is still 
not proportional to the third one which shows that span e (A® 2 ) ^ core® (A). However, it 
can be checked that this happens in span ffi (v4® 4 ), with the first two columns still equal to 
the generators of V($(A® 2 , 1), which shows that span e (74® 4 ) is the sum of above mentioned 



max cones, and hence span ffi (y4 $ 



spam 



core® (A). Hence we see that 



A is column periodic (S§) and the core finitely stabilizes. See Figure |2] for a symbolic 
illustration. 

In nonnegative algebra, p + = 1 and p+ = 0.7924. Computing the eigenvectors of A 
and A x2 yields 

V+(A, 1) = span + {(0.1326, 0.1326, 0.6218, 0.7604)}, 



V+(A 



x2 



span + {(0.2646, 0, 0.5815, 0.7693), (0, 0.2566, 0.6391, 0.7251)}, 



and 



V+(A,p+) = span + {(0, 0, 0.6612, 0.7502)}. 
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Figure 2. The spans of matrix powers (upper graph) and the periodic 
sequence of their eigencones (lower graph) in Example 2 (max algebra) 

Here core + (A) must be equal to the ordinary (Minkowski) sum of V + (A x2 , 1) and 
V+(A, Py). To this end, it can be observed that, within the first 4 digits, the first two 
columns of A xt become approximately periodic after t = 50, and the columns of powers of 
the normalized submatrix A vv j p^ approximately stabilize after t = 7. Of course, there is 
no finite stabilization of the core in this case. However, the structure of the nonnegative 
core is similar to the max-algebraic counterpart described above. 
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